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The linear term proportional to \N — Z\ in the nuclear symmetry energy (Wigner energy) is 
obtained in the framework of a model that uses isovector pairing on single particle levels from a 
deformed potential combined with a interaction. The pairing correlations are calculated by 
numerical diagonalization of the pairing Hamiltoriian acting on several levels nearest the Fermi 
surface. The experimental binding energies of nuclei with N Z are well reproduced without 
introducing new parameters. 
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The nuclear ground state energy, E{N, Z), as a func- 
tion of the proton number (Z) and neutron number (N) 
or atomic mass number {A = N + Z) is very well de- 
scribed by the celebrated empirical mass formula (see 

e.g. n) 

E{N,Z) = Ev+Es+Ec+Ea+Ew+Ep+Eshell- (1) 

The various terms have a clear physical meaning. The 
volume term, Ey = —ay A, describes the constant bind- 
ing energy of a nucleon in saturated nuclear matter. The 
surface energy, Eg = agA^^'^, accounts for the lack of 
neighbors in the surface. The term, Ec — acZ"^ j A^l'^ ^ 
is the electrostatic Coulomb energy. The (a-)symmetry 
energy, E^ — a^(iV — consists of two approxi- 

mately equal contributions. The "kinetic" part accounts 
for the Pauli principle, which requires the nucleons to oc- 
cupy higher single particle levels with increasing asym- 
metry \N — Z\. The "interaction" part originates from 
the difference between the proton-proton and neutron- 
neutron interactions and the proton-neutron interaction. 
The pairing energy (Ep) describes the energy gain by 
forming pairs of protons or neutrons. The shell energy 
{Eshell) is a manifestation of the level bunching around 
the Fermi level. The term, E^r — aw\N — Z\/A, is called 
Wigner energy, because Wigner [2] gave a first interpre- 
tation in terms of his super multiplet theory. However 
its physical origin has been the subject of a long debate, 
which has been recently reviewed by [3]. Modern mean 
field approaches reproduce the ground state energies very 
well, except the Wigner energy, which has to be added 
as an ad- hoc phenomenological term (see e. g. This 
means that the physics behind the Wigner energy is not 
taken into account by present mean field theories. In 
this letter we demonstrate that the Wigner energy is 
obtained, without introducing any new parameters, by 
including the isovector proton-neutron pair correlations 
determined by numerical diagonalization of an isorota- 
tional invariant pairing Hamiltonian. 

Experimentally, the coefficients aA and aw are not 
very different. As the ground state isospin (T) of most 
nuclei is equal to their isospin projection (T^ = =^^-5^). 



The sum of the symmetry and Wigner energies is ap- 
proximately proportional to T(T -\- 1). The T - depen- 
dence is suggestive, because the isospin operators obey 
the same SU2 algebra as the angular momentum opera- 
tors. Spontaneous breaking of the rotational symmetry 
by the deformed mean field leads to the appearance of ro- 
tational bands. The energies of the rotational levels are 
proportional to /(/-l-l), with / being the angular momen- 
tum. The analogy between nuclear spin and isospin led 
Fraucndorf and Sheikh to suggest that the T(T+ 1) 
dependence of the ground state energy is a manifestation 
of an isorotational band. 

The band appears because the isovector pair field, 
which is a vector, spontaneously breaks rotational sym- 
metry in isospace. Satula and Wyss discussed the anal- 
ogy of the cranking model in isospace [71 [8] . In the limit 
of strong symmetry breaking, simply the isorotational 
energy T{T + l)/20 is added to the intrinsic energy of 
the symmetry breaking mean field, the orientation of 
which can be taken such that the proton-neutron pair 
field is zero [H |6j. Refs. [9| and |in| used this simple 
limit for successfully interpreting the excitation spectra 
of nuclei with N ^ Z. Applying the Mean Field and Ran- 
dom Phase Approximation to an isorotational invariant 
isovector pairing interaction, Neergard has reproduced 
the experimental observation a a ~ aw [21 [HI [l2] . The 
virtue of such an approach is that the Wigner energy ap- 
pears without introducing any new parameter, because 
the strength of the proton-neutron pair correlation is 
fixed by the isorotational invariance of the isovector pair- 
ing Hamiltonian. Both approaches only work well when 
the isorotational symmetry is substantially broken by the 
isovector pair field. However, this not always the case of 
the medium mass nuclides, where the Wigner energy is 
seen to vary. For this reason, we take the pair correlation 
into account by numerically diagonalizing the isovector 
pair Hamiltonian within a configuration space spanned 
by seven single particle levels nearest the Fermi surface. 

As a starting point we assume that the isospin mix- 
ing caused by the Coulomb interaction can be neglected. 
Ref. |13j estimates the admixture of components with 
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T > Tz to the ground state to be of the order of 0.3 % 
for A ~ 70. With this assumption, the Coulomb energy 
can be separated from the energy caused by the strong in- 
teraction. Following [13] (and earlier work cited therein) 
we subtract the Coulomb energy from the experimen- 
tal energies and compare the resulting energies with our 
model. The Myers-Swiatecki expression for Ec [IS] was 
used to fit the Coulomb energy difference AEq between 
the experimental binding energies [TB] of 75 pairs of even- 
A mirror nuclei mirror nuclei for 20 < A < 100. We find 
that: 
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= 0.6967^2/3 _ 0.8027[MeV], 



(2) 



which is comparable to previous fits (see e.g. [TJ). After 
generalizing the restriction for mirror nuclei, the resulting 
expression for the Coulomb energy is determined: 

Ec = 0.697-— - 0.240— - 1.064— -[MeF]. (3) 
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In the following we use "experimental strong interaction 
energies" defined by Es = Eexp{N, Z) - Ec{N, Z). 

In accordance with the concept of isorotational 
bands, we write the energy of an isobaric chain 
{A = Z ^ N ~ const) in the form 



E{N, Z) = E,, 



T{T + X) 
26 



r = IT, 
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where, Eint is the energy of the intrinsic {N = Z) config- 
uration. As discussed in [51 [3], the ordinary BCS ground 
state without proton-neutron pairs is a legitimate intrin- 
sic state. It is a mixture of only even N and Z, which 
implies that must be even if A/2 is even or must be 
odd if A/2 is odd. Hence the ground state isorotational 
bands of even-even nuclei are composed of even values of 
T = if A/2 is even and odd values of T = if A/2 is 
odd. The term 1/2G is a combination of the coefficient 
as of the the symmetry energy and a contribution from 
the shell energy Eshell{Z,N), which depends on T^. 
Likewise, X/Q is a combination related to the coefficient 
aw of the Wigner energy, which also contains a contri- 
bution from EsHELL- We introduce the experimental 
isorotational frequency 
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E{T + 2) - E{T) _T + 1+X 
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The slope and intercept with the w-axis determine 1/0 
and X. We take the experimental ground state energies 
of the three nuclei =0, 2, 4 if A/2 is even and of Tz =1, 
3, 5 if A/2 is odd and calculate two points of lo{T) be 
means of Eqn. [5] Note that this is just a recombination 
of the experimental ground state energies, which aims at 
exposing the the Wigner energy. Figs. [2] and [3] show the 
experimental values of 1/0 and X, where the experimen- 
tal binding energies are taken from the most recent mass 



evaluation from [TB]. The small error bars for the lower 
mass region are primarily caused by the uncertainty in 
the spherical Coulomb energy fit. The large error bars 
in the A > 80 region are mainly caused by the error in 
binding energy of the nucleus nearest to or at N = Z . 

A monopole isovector Hamiltonian is used to describe 
the pair correlated ground state. 
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where pj^ and create a proton and a neutron, respec- 
tively, on the level k, and k denotes the time reversed 
state of k. This Hamiltonian is invariant under rotations 
in isospace, i.e. it conserves isospin. The pairing problem 
is solved via matrix diagonalization for the three nuclei 
with T^=0, 2, 4 or T^^l, 3, 5. 

As the computational complexity grows quickly with 
number of levels, the diagonalization is carried out within 
the configuration space spanned by the set of single parti- 
cle energies tk centered around Fermi level of the N = Z 
nuclide within the considered isobaric chain. The matrix 
is constructed by successive application of the pairing in- 
teraction onto the uncorrelated ground state, which en- 
sures that the configuration space contains only states 
with T = Tz- Shell effects in the non-interacting lev- 
els outside of this window are accounted for by sum- 
ming the single particle energies for all occupied pro- 
ton and neutron levels. The single particle energies are 
calculated by means of the micro-macro method using 
a Nilsson Hamiltonian as described in Ref. [18 . For 
each nucleus the equilibrium deformation has been cal- 
culated. In these calculations BCS pairing was used with 
Az = liAMeV/A^/^ and Ajv = UmieV/A^'"^ as sug- 
gested in |19| , which means the deformation of the intrin- 
sic state is found. Ref. [23 discusses this procedure in 
more detail. The resulting deformations are comparable 
with those from Ref. [21]. The single particle energies 
used in the diagonalization of the pairing Hamiltonian are 
taken as the average of the proton and neutron energies 
calculated by the Nilsson model at equilibrium deforma- 
tion. 

The pairing strength G must be adjusted to match the 
reduced single particle space. We use experimental stag- 
gering 2Aeeoo{N, Z) between the lowest T = states 
in the even-even and odd-odd N = Z nuclei to determine 
GiA). 



2Aeeoo{N,Z) = 
\E(N - 1, Z - 1) - 2E{N, Z) - E{N + 1, Z + 1)| 



(10) 



3 




-•t ■ 1 ■ 1 1 1 1 1 ' 1 1 1 ■ 1 ■ r O-^i . 1 . 1 . 1 1 1 1 1 . 1 . 1 . r 

20 30 40 50 60 70 80 90 100 20 30 40 50 60 70 80 90 100 

Atomic Mass Number (A) Atomic IVIass Number (A) 



FIG. 1: The experimental even-even to odd-odd energy differ- 
ence 2^EEOO compared with a seven level pairing calculation 
with the interaction strength (111. 



The same quantity was calculated from the theoretical 
energies obtained by diagonalization. The lowest T = 
state in the odd-odd nucleus has two-quasiparticle char- 
acter (see [5]), i. e. its Fermi level is blocked. Accord- 
ingly, we remove it from the diagonalization (as discussed 
in Pl]), which is carried out for iV = Z = A/2 - 1. Two 
times the blocked single particle energy is added to the 
resulting energy. Fig. [l] shows that the calculations re- 
produce the experiment rather well using the following 
pairing strength: 

7 fi'^ 

G=^,\MeV\. (11) 

The use of a fixed set of single particle levels within 
an isobaric chain only accounts for the "kinetic" part of 
the symmetry energy, which results in underestimating 
it by roughly a factor of two. Neergard suggested adding 
the schematic interaction CT ■ T in order to obtain the 
correct symmetry energv[ll). On the mean field level, 
this "symmetry interaction" generates a shift of the sin- 
gle particle energies by C(N — Z)tz, which mimics the 
realistic isospin dependence of the single particle ener- 
gies responsible for the interaction part of the symmetry 
energy (see [13] for a more detailed discussion) . Since 
T is a good quantum number, the symmetry interaction 
only adds CT(T + 1) to the calculated energies. Fig. [2] 
demonstrates that 

24 Q4 

C^—[MeV] (12) 

provides an excellent reproduction of the experimental 
values of 1/6. The fluctuations are caused by the irreg- 
ular spacing of the single particle levels, i. e. the shell 
structure. They are more pronounced for the spherical 
nuclei, which have strongly bunched levels, whereas for 
62 < A < 90 nuclear deformation generates more evenly 
distributed single particle levels. 



FIG. 2: The curvature of the symmetry energy (1/^) com- 
pared to experiment. 

Fig. |3] demonstrates that the calculated Wigner X well 
reproduce the experiment. As in the case of 1 /0, the fluc- 
tuations reflect the irregular level spacings. Varying the 
coupling constant G, we found that X approaches the 
strong symmetry breaking limit of 1 only for G > 2MeV 
in the case of spherical nuclei. In the case of equidis- 
tant level it is reached for G > 0.5MeV. Eq.(ll) gives 
O.SMeV for A = 80, i. e. the pairing correlations are not 
strong enough to average out effects of the substantial 
level bunching. For most of the isobaric chains studied, 
the calculated X shows the same up and down varia- 
tion, which is seen experimentally. The size of fluctu- 
ations of X turned out to be sensitive to the level po- 
sitions near the Fermi level, which cannot be expected 
to reproduced in detail by our mean field. The calcula- 
tion yields X values which are too small in the region 
72< A < 92. A possible reason for the discrepancy may 
be the theoretically determined deformations, which are 
fairly constant and small, whereas the experimental in- 
formation (_B(£'2,2^ — >■ 0^), rotational spectra) points 
to large deformations for some of the considered nuclei. 
A supplemental calculation intended to simulate changes 
in deformation along an isobaric chain, as indicated by 
the experiment, results in X values around 2. 

Odd-odd N — Z nuclei with A > 40 often have a 
ground state that has T = 1 > = 0. The inversion of 
the isospin order has been explained by Refs. [23] and 
[24] . The T = state in the odd-odd nucleus is lifted rela- 
tive to the r = ground state of the even-even neighbors, 
by the two-quasiparticle excitation energy 2Aeeoo- The 
T = 1 state is lifted by the isorotational energy 1/0, 
which is somewhat smaller than 2Aeeoo- Fig. ]4] shows 
that calculations fairly well reproduce energy difference 
between the lowest T — 1 and T — states, which mea- 
sures to relative strength of the pair correlation and the 
isorotaional energies. 
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FIG. 3: The Wigner X compared to experiment. 
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FIG. 4: Energy difference of the first T = and first T = 1 
states in odd-odd N=Z nuclei calculation compared to |25j 
experiment. 

We have demonstrated that a model based on single 
particle levels in a deformed potential, isospin conserv- 
ing isovector pairing, and a schematic " symmetry" inter- 
action proportional to reproduces the term linear in 
|iV — Z| (Wigner energy) in the nuclear binding energy. 
The pairing correlations are treated exactly, not in mean 
field approximation, by numerical diagonalization in a 
space of seven single particle levels. The model does not 
introduce new parameters as compared to standard mean 
field approaches: The pairing strength is fixed by the 
even-even to odd-odd mass difference, and the strength 
of the symmetry interaction is determined by the symme- 
try energy. The shell structure causes strong fiuctuations 
of the Wigner energy which are fairly well described by 
the model. The remaining deviations can be attributed 
to inaccuracies of the single particle energies. The results 



seem to provide no evidence for the presence of isoscalar 
proton-neutron pair correlations, in contrast with claims 
in Refs. [HIH]. Our results suggest that a combination 
of an isorotational invariant effective interaction in the 
particle-hole channel with isovector pairing interaction 
will give the Wigner energy, provided the pairing corre- 
lations are treated beyond the mean field approximation 
and isospin is conserved. How to accomplish this for the 
present standard mean field approaches remains to be 
studied. In a future study we will address this question 
by comparing our results with approximations as e. g. 
isospin projected mean field solutions. 
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